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1 Introduction 



Field-theoretical calculations in the context of critical phenomena [1,2] fall 
into two categories (i) those based on expansions about the upper [3] or lower 
[4] critical dimension (e-expansions) and (ii) those that are carried out at fixed 
dimension [5]. Within these two approaches, one can further distinguish be- 
tween two types of renormalization: the use of renormalization conditions [6] 
and the minimal subtraction scheme [7]. Whereas the e-expansion has been 
used extensively with both types of renormalization, the theory at fixed dimen- 
sion is usually presented within the framework of renormalization conditions 
[1,5,8-13]. 

In this paper we shall use an approach [14-17] that combines the theory at 
fixed dimension with the minimal subtraction scheme. While any renormaliza- 
tion scheme implies a (scheme- dependent) decomposition of correlation func- 
tions into exponential parts and amplitude functions our present approach is 
especially advantageous since here the exponential parts are determined en- 
tirely from pure dimensional poles which remain unchanged in extensions of 
the theory from T > T c to T < T c as well as to finite k, u and L (wavenum- 
ber, frequency and system size) [14-21]. Accordingly, this concept of a min- 
imally renormalized theory at fixed dimension has been applied successfully 
not only to static critical phenomena [22-24] but also to dynamics in equi- 
librium [15,20,25-27] and nonequilibrium [28], finite-size effects [29-31] and 
surface critical phenomena [32,33]. Here, we use this approach to calculate 
static amplitude functions of 0(n) symmetric systems below Tc- For n > 1 
these systems are of particular interest because of the massless Goldstone 
modes [34] governing the long-distance properties. 

The 0(n) symmetric 4 model [1,2] describes the most common examples of 
critical behavior: in liquid-gas systems (n — 1), superfluid 4 He (n = 2) and 
isotropic magnets (n = 3). Since the critical properties of 4 He near the A- 
transition can be measured with high accuracy not only near a single critical 
point T\ [35] but also near a whole (pressure-dependent) line T\(P) of crit- 
ical points [36] the opportunity exists for the universality predictions of the 
renormalization group (RG) theory to be tested at a highly quantitative level. 
Owing to the wide range of applicability of the RG concept [1] this test is of 
relevance not only to statistical physics but equally well to elementary par- 
ticle physics and condensed matter physics. In an effort to obtain the most 
accurate data possible, high-precision measurements are to be carried out in a 
microgravity environment [37]. The corresponding theoretical challenge is to 
calculate as accurately as possible the properties of the 4 model appropriate 
for a comparison with the helium data. This includes not only the well-known 
exponent functions [1,2,16] whose fixed points values determine the critical ex- 
ponents but also the less well-known amplitude functions [14-17,22,23] which 
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determine the equation of state and which contain information about asymp- 
totic ratios of leading and subleading amplitudes [38]. They are also needed 
for a description of nonasymptotic critical properties. The latter are impor- 
tant in order to distinguish universal from non-universal properties of critical 
behavior over a wide temperature range [19,39]. 

Accurate RG calculations of amplitude functions for both T > T c and T < T c 
have been performed within the d = 3 theory only for the case n — 1 [10- 
13,22,23]. In this case Borel resummations yielded accurate results because 
the perturbation series are known to high enough order (five loops). For n > 
1, only amplitude functions above Tc have been obtained with comparable 
accuracy [10,22]; below Tc, the complications due to Goldstone singularities 
have prevented calculations to sufficiently high order for resummation methods 
to yield accurate results. For example, the amplitude functions for the (n > 1) 
order parameter and superfluid density have been computed in the d = 3 
theory only in one-loop order [17,39] where Goldstone singularities do not yet 
arise. 

Our aime here is to present the first step towards filling this gap of theo- 
retical knowledge for n > 1 below T c for the O(n) symmetric 4 model in 
three dimensions. Specifically, we calculate in two-loop order the amplitude 
functions of the order parameter, the specific heat and the helicity modulus 
introduced by Fisher, Barber and Jasnow [40] . We also calculate the q 2 part of 
X T (?) _1 where x T (q) is the wavenumber dependent transverse susceptibility. 
This quantity enters Josephson's definition [41] of the superfluid density and 
provides an independent check of our calculation of the helicity modulus. 

These calculations serve three purposes. First, by comparing the zero-, one- 
and two-loop terms we can study the reliability of the low-order perturbation 
theory. Whereas the two-loop contribution to the q 2 part of x T (g) _1 turns 
out to be very small we find that the two-loop contributions to the superfluid 
density and specific heat below T c are comparable in magnitude to the one- 
loop contributions, thus indicating the necessity of higer-order calculations 
and Pade-Borel type resummations. Second, we can anticipate some of the 
technical difficulties that will also appear in higher-order calculations, where 
the amount of computational labor required is substantially greater. These 
difficulties include the removal of ultraviolet divergences in three dimensions 
and of spurious infrared (Goldstone) divergences, both of which first appear 
here at the two-loop level. Third, we provide part of the information for a 
future study of the n-dependence of amplitude functions below T c beyond one- 
loop order to be carried out in a combined analysis of known Borel-summed 
amplitude functions for n — 1 and of low-order contributions for n > 1. 

This paper is organized as follows. In Sec. 2, we summarize some relevant 
aspects of the theory at fixed dimension. In Sec. 3, because of our interest in 
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computing the helicity modulus (superfluid density), we study the free energy 
for a state with a twisted order parameter. In Sec. 4, we use the free energy 
to derive the bare two-loop expression for the equation of state and the lon- 
gitudinal susceptibility. We follow this, in Sees. 5 and 6, by a calculation of 
the specific heat, of the helicity modulus and of the q 2 part of x T (g) _1 . We 
find several contributions (diagrams) that have Goldstone divergencies which 
we regulate by use of an external field. For rotationally invariant quantities 
we show that these divergencies cancel leaving a finite result in the limit of 
zero field. For the order parameter, specific heat, helicity modulus (superfluid 
density) and for the q 2 part of x T (q , ) _1 , we determine their respective ampli- 
tude functions up to two-loop order. The results and conclusions are presented 
in Sec. 7 and details of the calculations are given in the appendices. Further 
details are given in Ref. [42]. A short summary of some of the results of the 
present paper has been presented in Ref. [43]. 



2 Field theory in three dimensions 



In this section, we outline the strategy of our calculations within the minimally 
renormalized massive field theory at fixed 2 < d < 4. Consider the 0(n) 
symmetric 4 model in the presence of an external field ho(x), as described 
by the standard Landau-Ginzburg- Wilson functional [1,2] 

H{0 o (x)} = J d d x Qr„$ + \ E(V0o,) 2 + «o(0o) 2 - ho ■ $o) , (1) 



where 0o( x ) — (0oi( x ), • • • >0On(x)) is an n component vector subject to the 
statistical weight ~ exp(— H) and V is the volume. The spatial variations of 
0o ( x ) are restricted to wavenumbers less than some cutoff A. The Gibbs free 
energy per unit volume (divided by ksT) is 

F (r , u , {Ko(x)}) = -V~ l In f D0 o (x) exp(-H{0 o (x)}) (2) 



and is related to the Helmholtz free energy T per unit volume via 



r (r , uo, {(0o(x))}) = F (r , u , {K (x)» + V' 1 J d d x £ (x) • (0 o (x)> , (3) 

v 

(4) 



<<Mx)>= -v -=<-?- 

dh (x) 
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Tq is the generating functional for vertex functions and can be obtained per- 
turbatively from the negative sum of all one-particle irreducible (1PI) vacuum 
diagrams. We shall consider always the bulk limit V — > oo. We shall fur- 
ther suppose that finite-cutoff effects are negligible (although this may not be 
justified in certain cases [21,44,45]) and that all integrals are evaluated in the 
limit A — > oo according to the prescriptions of dimensional regularization [1,2]. 
The ultraviolet divergences are thus manifested as poles in the dimension d 
at discrete values d\ = 4 — 2/1 (where / = 2, 3, . . .) if the vertex functions are 
considered as functions of r [10,11,16,46]. Our ultimate goal is to calculate 
the vertex functions at d = 3 as functions of the correlation length as this will 
ensure that the perturbative expansions have no poles and are (presumably) 
Borel resummable [10,11,16,17]. Our strategy will be to treat the d — 3 poles 
at the level of the bare free energy and use the resulting finite expression to 
derive bare expressions for the other quantities directly in three dimensions. 

One way to remove the poles at di < 4 is to rewrite the perturbation series in 
terms of the variable tq — tq c where ro c is the critical value of ro = tq c + a^t 
with t — (T — T c )/T c being the reduced temperature. The structure of r 0c 
is [46] r 0c = u^/ e S(e), where S(e) is a dimensionless function which is finite 
for e = 4 — d > except for the poles at d — di (e = 2/1) which cancel 
the poles mentioned above. The role played by r 0c in the d = 3 theory was 
discussed in detail by Bagnuls and Bervillier [10], by Bagnuls et al. [11] and, 
in the context of the minimal renormalization, by Schloms and Dohm [16]. 
Another possibility, which we shall use below, is to use a shifted variable r' 
(given in Eq. (11) below) which differs from r — r 0c only by a conveniently 
chosen finite constant [11,23]. Expressed in terms of the variables r — r 0c or 
To, the unrenormalized vertex functions are finite for d < 4 but still have poles 
at d = 4 which can be absorbed by the Z-factors either by use of (d = 3) 
renormalization condititions [10,11] or within the minimal subtraction scheme 
[14-16]. Although in the latter approach the Z factors are determined by 
(d = 4) pole terms ~ e~ n this does not imply the necessity of an e expansion, 
as shown in Ref. [16]. The Z-f actors connect the bare and the renormalized 
quantities in the usual way 

r = Z r _1 (r - r 0c ), u = iT € A d Z~ x Z\u^ $= Z^ 1/2 0o, (5) 
where in the minimal subtraction scheme to two-loop order [1,2] 



16(n + 5)^r -20- 
e e 

> + 8) 2 - - (5n + 22)- 
e z e 



Z r (u, e) = 1 + 4(n + 2)-u + (n + 2) 
Z u (u, e) = 1 + 4(n + 8)-u + 16 
Z (w,e) = l-4(n + 2)-U 2 + 0(u 3 ). (8) 



u 2 + 0(u 3 ), (6) 
u 2 + 0(u 3 ),(7) 
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Fig. 1. Two- loop diagram of the two-point vertex function T Q ' above T c whose 
dimensionally regularized integral has a pole at d = 3 as contained in Eqs. (A.l), 
(A. 5) and (A. 7) of App. A. This pole is also contained in similar diagrams of T 
and r y rf in Figs. 3 and 4. 

In Eq. (5), fi^ 1 is a reference length and 

A d = r(l + e/2)r(l-e/2)S' d (9) 

is a convenient geometrical factor where T is Euler's gamma function. For 
applications to amplitude functions at d < 4, this factor Ad has considerable 
advantage [14-17] over the more commonly used [1,2] geometrical factor Sj = 

[2 d - 1 7r d / 2 r(rf/2)]- 1 . 

The pole at d — 3 in r 0c can be traced to two-loop diagrams of the type shown 
in Fig. 1. Having determined its coefficient (see Appendix A), we write 

r 0c = - 2 {n + 2)^- + S(e, n) uf (10) 

where S(e,n) is finite for e = 1. As shown previously [16] an explicit knowl- 
edge of the function S(e, n) is not necessary since S(e, n) does not enter the 
final expressions of thermodynamic quantities as functions of the reduced tem- 
perature t. For convenience, therefore, we use instead of r — r 0c the variable 
r 'o = r o — $ r o where 5r contains the d — 3 pole of r 0c but not the poles of 
S(e, n) at d ^ 3. Thus we take 

r'o = r - ~ 2 {n + 2)^1 _ C (n) uf (11) 

with a conveniently chosen finite constant C{n) whose value is given in Eq. (A. 9). 
For n — 1, this r' coincides with the variable r' used previously [11,23]. The 
final results do not depend on the choice of C(n). 

When the perturbation series are expressed in terms of the variables r — r 0c 
(or To ) the series at d — 3 contain logarithms of u arising from the nonan- 
alytic -uo dependence of tq c (or Stq). For higher-order calculations intended 
for resummation, such logarithms are inconvenient since the series are then 
not Borel-summable. Perturbation series that are free of these logarithmic 
terms (and which are also finite in d — 3) can be obtained if they are instead 
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expressed in terms of the correlation length [16,17] (or some other physical 
quantity). 



Above T c , the correlation length £ + is defined as usual by 



X+(0)9x + (?)- 1 /^ 5 



q=0 



;i2) 



where x+(q) is the bare susceptibility at finite wavenumber q. Below T c , a com- 
mon definition of the correlation length for n — 1 and n > 1 is complicated by 
the fact that, for n > 1, the spatial decay of the order parameter correlations 
is not exponential [38,40,47,48]. Therefore, it is not straightforward to define 
an analogous quantity which plays the same role as above T c in removing 
the logarithms in uq from the perturbation series obtained below T c . Here, 
we follow Schloms and Dohm [17] who have introduced a pseudo-correlation 
length £_ which performs precisely this task. We use Eq. (12) and the defini- 
tion of given in Ref. [17] to determine the two-loop expressions for 
function of £ + and £_ in three dimensions. The results are (see Appendix A) 



r' = Z+ 2 { l + (n + 2) 



-u £ + + ^(u £ + ) 2 ( — + 21n(24 Mo £+) 



(13) 



-2r' = T 2 | 1 + (n + 2) [i^- - ^«-) 2 + 41n(24^- 

+o(^e)j 



(14) 



for T > T c and T < T c , respectively. For n = 1, Eqs. (13) and (14) agree 
with Eqs. (All) and (3.8) of Halfkann and Dohm [23]. These formulas are not 
to be regarded as relations between the correlation length and the reduced 
temperature t but rather as an intermediate step in the calculation of the 
bare quantities as functions of the correlation lengths £±. The bare quantities 
can finally be obtained as functions of t via the connection between £± and t 
as given in Appendix A. 

For the purpose of deriving physical quantities as a function of £± in three 
dimensions it is not nessecary to calculate the perturbation series first as a 
function of r f at d — 3 and then to substitute r' as a function of £±. A more 
direct way would be to start from (unevaluated) diagrammatic expressions as 
a function of r at d ^ 3 and to substitute r in the form of (unevaluated) 
diagrammatic expressions as a function of £± at d ^ 3. The resulting perturba- 
tion series would consist of integral expressions that have both finite limits for 
d — > 3 and are free of logarithms in uq. The advantage of this procedure is that 
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it avoids an explicit treatment of d = 3 poles and it requires an evaluation only 
of a simplified form of finite integral expressions at d — 3. We have not chosen 
this more direct route of calculation in the present paper since at the two-loop 
level its advantage is not yet substantial. In future calculations beyond two- 
loop order, however, the simplification and reduction of computational labor 
implied by this procedure may become crucial. For present purposes, the in- 
termediate expressions given in terms of r' illustrate the appearence of the 
lnu terms and enable us to make contact with the earlier work of Refs. [11] 
and [23]. 

Finally, let us note that there is some flexibility in the definition of for 
general n. While it would be natural to define so as to coincide, for n — 1, 
with the usual correlation length of Ising-like systems below T c , a suitably 
modified definition of for n > 1 could well absorb the poles at d < 4 
and the logarithms in u and yet incorporate higher-order terms that are 
better adapted to the region T < T c than those of of our paper. (The 
higher-order terms of our £_ in Eq. (14) are determined essentially from x+(q) 
above T c , see Appendix A.) A definition of £_ for n — 1 via the q 2 part of 
X-iq)" 1 could formally be extended to n > 1, for example, by including only 
those diagrammatic contributions that cause an exponential spatial decay of 
the correlation function (i.e., by omitting those transverse parts of diagrams 
causing the algebraic decay). This may lead to a simplified representation of 
the amplitude functions of other physical quantities below T c at higher order. 



3 Bare free energy 

Because of our interest in computing the superfluid density for 4 He (n = 
2), we consider a state in which the order parameter (0 O ( X )) has a uniform 
twist [40,49] along a fixed direction specified by a wave-vector k. For n — 2, 
this twisted state is equivalent to the situation in 4 He at constant superfluid 
velocity v s = Kk/m where the order parameter is a complex macroscopic 
wavefunction of plane wave structure (0o( x )) — Vo( r o, exp (ik-x). In three 
dimensions this case has been studied recently in one-loop order [24]. Here, 
however, the situation is more delicate than in the earlier work on account 
of the (spurious) Goldstone divergences that plague the perturbation theory 
beyond one-loop order. 

We begin by introducing, for n > 2, an external field 

— * 

ho(x) = (/i cosk-x, /i sink-x, 0, ... ,0) (15) 
which not only regulates the Goldstone divergencies but also generates the 
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twist in the order parameter; the amplitude ho is independent of x. At the 
same time, we introduce a local rotation according to 



0oi(x) = </>oi( x ) cosk-x — (x) sink-x, 
002 (x) = O1 (x) sink-x + O2 (x) cosk-x, 

Ol (x) = Ol (x), 7^ 1,2). (16) 

Substituting Eqs. (15) and (16) into Eq. (1), we obtain 



H{0 o (x)} = H'{0 o (x),k} = J d d x 



ir o 0o 2 + ^ 2 (0 o 2 1 + o 2 2) + ^E(V0 Ol ) 2 



+ k.j + Mo (0o 2 ) 2 -^o-4 



(17) 



where j(x) = ol (x)V0 O2 (x) — O2 (x)V0 ol (x). The advantage of working with 
the variables O j(x) is that o (x) is coupled to the spatially homogeneous field 
h'o = (h , 0, . . . , 0)' through h' • O = h (f)' 01 . We anticipate that the order 
parameter (0 O ) = (M , 0, . . . , 0)' in the rotated system will also be homoge- 
neous. In the original coordinate system, this corresponds to the twisted order 
parameter [40,49] 

(0 o (x)) = (M cosk-x, M sink-x, 0, . . . , 0) (18) 



where M is independent of x; the k-dependent terms in Eq. (17) represent the 
additional energy associated with the applied twist. This is the same situation 
as studied previously [24,40,49] except that now M (r , k, h ) depends not only 
on r and k but also on ho. Since the twist affects only the orientation of the 
order parameter in the 1-2 plane of its n-dimensional space, there are (n — 2) 
equivalent transverse components left. Their equivalence will be manifested 
through terms proportional to (n — 2) in the expression for the free energy in 
Eq. (26) and Fig. 3. 

We are now in a position to set up the perturbation theory. We make an 
expansion around the exact average (0 O ) according to 

0o(x) = (0o)+50 o (x) (19) 



— * 

so that (50o (x)) = 0. Substitution of Eq. (19) into Eq. (17) and use of the 
Fourier representation 

50 o (x) = J ^(P)e ipx (20) 
p 
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where J p = (27r)' d J d d p, leads to 



H'{0 o (x),k} = H'{(0 o ),k} + J d d x 



50 o (x) 



1 



+ -J ^ (-p)K(p,k)Vo(p) 
p 

+ I d d x [4 Mo ((0 o )-50 o ) ( 50 o 2 + % (50 o 2 ) 2 -/; o -50 o 
where K is the n x n matrix 

/*T(p,k) 



(21) 



tf(p,k) = 



V (f T+p 2 )/n-2 



(22) 



J„__2 is the (n — 2) x (n — 2) unit matrix and 



tf(p,k) 



r L + & 2 +£> 2 2ikp 
— 2ikp f T + A; 2 +p 2 



r L = r + 12m M 2 , r 0T = r + 4m M 2 



(23) 



(24) 



(see also Ref. [24]). For k = 0, the diagonal elements of K(p,0) 1 yield the 
standard longitudinal and transverse propagators 



G L (p) = (r 0L + PT , G T {p) = (f 0T + P 



2\-l 



(25) 



The diagrammatic rules corresponding to Eqs. (21)-(25) are indicated in 
Fig. 2. The difference between the present perturbation theory and that for 
k = case lies in the k-dependence of the propagators and in the greater num- 
ber of vertices due to the three (rather than two) types of components for the 
order parameter. The resulting two-loop expression for the (bulk) Helmholtz 
free energy per unit volume reads (up to an unimportant additive constant) 



r (r , u , M , k) = - (r + A; 2 )M 2 + u M* + - J In [det K(p, k)] 

p 

+ i(n-2) J\n(f 0T +p 2 )+X (r ,u ,M ,k) + O(ul) (26) 
p 
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x =R L (p) - x-M^ = Rlt(p) 1 =G T {p)S aP 



2 2 
I X I 



R T (p) 



2 1 
— MX 



Rtl(p) 



1i 1 li 



1 / \ 1 2 / \ 2 «/ \ (3 

-4u M -4u M -4u M S a p 



L \ / 1 1 \ / 1 



1 / \ 1 2 / \ 2 «/ \ (3 

2 \ / 2 2 \ / 2 a\ / [3 



2 / \ 2 q/\/3 7/ \ J 

—mo — 2mo5q,/3 —uaSaPfS 

Fig. 2. Propagators and vertices corresponding to Eqs. (21)-(25). 5 a p is the Kro- 
necker delta and a, (3,^,5 = 3,4, ...,n. For i?^, Rt, Rlt, Rtl and see 
Eqs. (C.7), (C.8) and (25), respectively. 

where X is considered to be of 0{uq) and denotes the negative sum of the 1PI 
two- loop vacuum diagrams shown in Fig. 3. We label these diagrams A, B, . . . , 
L and refer to the associated integral expressions as X a, X ob , . . . , X 0L [see 
App. C]. We have omitted diagrams with tadpole insertions since their sum, 
being equal to (<50 o (x)), vanishes according to Eq. (19). Eq. (26) provides the 
basis for deriving the helicity modulus in Sect. 6. For k = 0, Eq. (26) reduces 
to the two-loop expression given by Chang and Houghton [50] and by Shpot 
[51]. 

Like the integral associated with Fig. 1, the integrals Xqh, Xqj and Xqj are 
ultraviolet divergent in three dimensions and thus have poles at e = 1. We 
denote these integrals by Xh, Xj and Xj when they are expressed as functions 
of r' , Eq. (11). In order to collect all d — 3 pole terms of r (r , u , M , 0) we 
consider X H , Xj and Xj at k = and define 



ulM%f(r' , u , M ) =lim 



^M^5r + X H + Xj + Xj 



(27) 

k=0 
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loop contributions reduce to products of standard one-loop integrals. Thus, we 
obtain for the bare Helmholtz free energy per unit volume in two loop order 
at k = and d = 3 



rK, I*,, M ) = l -r' Q M% + Mo M 4 - J-rjV 2 - J-( n - l)r % 2 

+ ulM 2 f(r' , M , M ) + 0(ul u 2 In u ) . (32) 

For the case n = 1, this expression reduces to the two-loop part of Eq. (3.18) 
of Bagnuls et al. [11] and of Eq. (3.3) of Halfkann and Dohm [23] and agrees 
with their coefficients F^k for 6 = 1 and 6 = 2. The latter agree also with 
those in Eq. (Al.l) of Guida and Zinn- Justin [13] for n — 1. Eq. (32) contains 
logarithms of the coupling via UqMq/^q, uo, M ) ~ O(uo\nuo) as expected. 
Unlike the case n — 1, however, the two terms proportional to n — 1 in Eq. (32) 
depend nonanalytically on r T through r^ 2 and r^ 2 . These nonanalyticities 
lead to perturbative terms in the derivatives of the free energy (with respect 
to Mo) that diverge when r T — > 0. The origin of these divergences are the 
well-known Goldstone modes [34] that result from the fact that, for a bulk 
homogeneous system with a rotationally symmetric order parameter, there 
is no restoring force to a global rotation of the order parameter at ho = 
[48]. For the rotationally invariant quantities that we consider here (square 
of the order parameter, specific heat, stiffness constant), the divergences of 
individual perturbative contributions should be spurious and should cancel 
among themselves leaving a finite result in the limit h — > 0, as required on 
general grounds [52,53]. Within our two-loop calculation, we shall show that 
this is indeed the case. This is in contrast to the transverse and longitudinal 
susceptibilities which possess true physical singularities in this limit [48,52,54- 
56]. 



4 Order parameter 



4-1 Bare theory 



The bare perturbative expression for the order parameter Mq at k = can 
be obtained either by evaluating the diagrams [51,57] of the one-point vertex 
function (see Appendix B) or, as we do here, by differentiation of the bare 
free energy T, Eq. (32), according to 

h (r' ,u ,M ) = JLt(r ,u ,M ). (33) 
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In three dimensions, we obtain 



1 3 

-rrho(/ , u ,M ) =r' + 4w M 2 

Mq 7T 



7T 



(n - l)«o»"Jr 



2tt 2 



(n - l)«o 



+ 2ulf{r' ,u ,M ) + —ul 



\r 0L J r 0L + 2(r 0T r 0L y/ 2 _ 

+ O(u 3 ,u 3 liiu ). (34) 



2 , io n r ° T 

n — n + 18 — 9 — 



— 1/2 

From the two-loop term ~ r 0T , it appears that this perturbative expression, 
in its present form, is problematic for small r T- In particular, Eq. (34) cannot 
be used directly at h = to obtain M . To circumvent this problem we invert 
Eq. (34) iteratively at h ^ 0. This leads to the bare perturbative form of 
the (implicit) equation of state (for Mg as a function of r' and h ) in three 
dimensions 



M, 



4m 



+ 



47T 



47T 



87T 2 



(n - l)u 



Qw 1 ' 2 + 3w(l + 2w 1/2 Y l -4m- 



1 Mn I q \ 1 ( , , , (-2r^, + 3x' 1 ) 1/2 
+ ^2-o(10 -n + 9w)- — 2 (n + 2)u In _ 

+ O(ul,ul\nu ) , (35) 

w = x -\-2r' + 3 X - 1 r 1 (36) 

where /i enters via the transverse susceptibility 

Xt = M /h (37) 



and w = O(h ) for h — > and r(, < 0. 

Eq. (34) also yields the bare longitudinal susceptibility Xz,( r o> u Oj M)) i n two- 
loop order since 



X, (r ,M ,M ) 



9 



<9M 



h (r' ,u ,M ) . 



(38) 



Combining this with Eq. (35), we find the perturbative form of x^ in three 
dimensions 



X: 1 = -2r' + 3*; 1 - ^o(-2r(, + 3X; 1 ) 1 / 2 [(n - 1)«T 1/2 (1 - 5w) - 3 - 9w 
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+ 



1 



47r 2 ~° 
4 



7T 



2 M 



n - ifw- 1 + 8(n - 1)«T 1/2 + (36 - 4ra - 5n 2 ) + W(iu) 

+ O(u 3 ,u 3 \nuo) (39) 



(n + 2) In — (n — 1) m 3 



where the function 



24w 



W(w) 



3m 1 / 2 



4(n - 1) + (8n - ll)w 1/2 - 6(n + l)iu - 30nw 3/2 

(40) 



(l + 2wV2)2 

-(30n + 42)u> 2 - 72u> 5/2 1 - 16(ra - 1) ln(l + 2u> 1/2 ) 



vanishes for w = 0. Note that the one- and two-loop terms in Eq. (39) diverge 
as w~ x l 2 and w~ l for ho — > below T c . This apparent breakdown of perturba- 
tion theory signals the Goldstone singularity of the longitudinal susceptibility 
[54-56] in this limit. An appropriate description of this singularity near T c is 
nontrivial [52] and will be studied elsewhere. 

The behavior of x L should be contrasted with that for the order parameter M . 
In deriving Eq. (35) from Eq. (34) or, equivalently, from the vertex function 



-1/2 



Tq (see App. B) we encounter terms below T c that diverge as x l J 2 ~ h ( 
for ho — > but which cancel among themselves. A similar cancellation of 
ln% T terms near four dimensions was noted previously by Brezin, Wallace 
and Wilson [58]. In three dimensions, we obtain from Eq. (35) the finite result 



o = ^(- 2 ^o) + ^-(-2<,) 1/2 - ^(n + 2)u In 



1 , _ f-2^) 1 / 2 



' /,: 8V ~' w ' 4tt v 2tt 2V " 1 24u 

+ -^H [10 - n + 4(n - 1) In 3] + O(u 2 , w 2 In w ) (41) 

for < and /io — > 0- The final step is to rewrite of Eq. (41) in terms of the 
correlation length £_ in order to remove the logarithms in uq. Using Eq. (14) 
for 7q as a function of we obtain the following two- loop expression for the 
square of the bare order parameter at h = 



M 2 (e_, Mo ,3) = r 1 



1 1 / 

+ — n + 8 - 



8w £- 8tt 
1 



-(1169ra+ 1042) 

864tt 2 v ; 



2tt 2 



n — 1) In 3 



(42) 



in three dimensions. For n = 1 this agrees with Eq. (3.13) and Table 2 of 
Halfkann and Dohm [23]. 
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4-2 Amplitude function of Mq in two-loop order 



The square of the renormalized order parameter M 2 is obtained from Mq (£_, u , d) 
according to [17] 



M 2 (£_, u, n, d) = Z^M^_, ^AfZ u Z^ 2 u, d) 

=e- d u(^,u,d) 



(43) 
(44) 



where is the amplitude function and where the renormalized parameters 
are defined in Eqs. (5)-(9). Solution of the RG equation for M 2 implies the 
following decomposition of the bare order parameter in three dimensions 



M 2 (£_, u , 3) = Z^u, 1) C 1 U(l, u(l-), 3) exp 



«(!-) 



[ UU ' ] du' 



(3 u (u, e) = -eu + (3{u) = u -e + fid^ In (Z^Zfy 



, (45) 

(46) 
(47) 



where and f3 u are the well known RG functions [16]. The flow parameter /_ 
is chosen according to = 1 and the effective coupling u{l) is determined 

by Eq. (A. 15). We recall that the exponential part in Eq. (45) is determined, 
within the minimal subtraction scheme, entirely by the d — 4 poles carried 
by Z^ and Z u and that the renormalization of Uq likewise involves only the 
pure poles contained within the product Z U Z^ 2 . Since these poles are the same 
above and below T c , we may use the Z-factors as given in Eqs. (6)-(8), without 
further calculation, to determine the amplitude function. In three dimensions, 
ft is given by 



U(l, u, 3) = Z^ 1 M 2 (£_, AnC'ZuZ^u, 3) 



(48) 



where we have used 



A 3 = (47T)- 1 



(49) 



and where Z^u, 1) and Z u (u, V)Z^{u, 1) 2 are given by Eqs. (7) and (8) eval- 
uated at e = 1. Using Eq. (42) we obtain from Eq. (48) 



3) = 



1 



327TU 



1 2 
— (160 - 82n) + -{n - 1) In 3 

277T 7T 



u + 0{u 2 ) (50) 



for the amplitude function in two-loop order. For n — 1, this agrees with Eq. 
(3.21) and Table 3 of Halfkann and Dohm [23]. Note that the vanishing of the 



16 



one-loop term in Eq. (50) is due to the choice of the geometric factor Ad in 
the definition of the renormalized coupling [17] in Eq. (5). 



5 Specific heat 



5.1 Bare theory 



Although the two-loop results for the amplitude functions of the specific heat 
(at ho = 0) have been given previously [14,17], we present here, for complete- 
ness, their derivation within the present approach. The bare specific heat C ± 
above (+) and below (— ) T c is given by 

C ± = C B - alt™ (51) 



where Cb is a noncritical background value and ao is a constant defined by 
(^o — r 0c) = dot. The vertex functions f ±'°^ can be obtained either from the 
sum of 1PI vacuum diagrams with two 0q insertions or by differentiation of 
the free energy, Eq. (32), according to 

fi 2,0) = ^jf(r , ,«o,M (r , o,tio)). (52) 



Using Eqs. (32) and (41) we obtain 



t™ = -^- 1/2 + 0(v>,v>\nuo), (53) 

f m = ~ ~ 7^" 2r or 1/2 + (" + 2)«o(-2r[ ) )- 1 + 0« u 2 lnt*>) 

(54) 

in three dimensions for r' Q > and r' < 0, respectively. As expected, Eqs. (53) 
and (54) contain no logarithms in u since, at this order, there are no diagrams 
of f ±'°^ with d = 3 poles. Nevertheless, in anticipation of future Borel sum- 
mations, we rewrite the above expressions for f j 2 ' '* in terms of the correlation 
lengths £± in Eqs. (13) and (14). Thus, we obtain 



f? 0) (£ + ,%,3W + 
fL 2 ' 0) (£-,«o,3)=£- 



16tt 32tt 2 
1 



-^- + A(n + 2)u ^ + O(uie) 



8w £_ 4vr 8tt 2 



(55) 
(56) 
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in three dimensions. In an earlier calculation [14] of the specific heat in three 
dimensions below T c (performed on the basis of the static distribution function 
of model C [59]) Goldstone singularities of individual diagrams were found to 
cancel among themselves, in accord with Eq. (54). 

5.2 Amplitude functions of C ± in two-loop order 

The renormalized vertex functions r^' - 1 are obtained from f±'°^(£±, u , d) ac- 
cording to [17] 

ri 2 ' 0) (£ ± , M ,M) = Z 2 r t^\i ± ,^Z u ZfAfu,d) - ^A d A(u,e) (57) 

where, in two-loop order, the RG function 

A(u,e) = -2n^ - 8n{n + 2)^u + 0{u 2 ) (58) 

absorbs the additive poles in both f^ 2 ' '' and f ( i' ' ) (see Eqs. (3.20) and (5.9) 
of Ref . [14] and Eq. (2.18) of Ref . [22]). Dimensionless amplitude functions F± 
can then be defined according to 

lt' 0) (£ ± , M ,M) = ~V- £ A d F ± (rt ± ,u,d) (59) 

so that the quantities of interest F±(l,u, 3) are given by 

F±(l,«,3) = -l^Zl^t^^A^ZuZ^u,?,) + A(u,l). (60) 

Use of Eqs. (55) and (56) together with the Z-factors Z r (u,l), Z u (u,l), 
Z<p(u, 1) and the additive function A(u, 1) of Eq. (58) leads to the two-loop 
formulas for the amplitude functions in three dimensions [14,17] 

F+(l, u, 3) = -n - 2n(n + 2)u + 0(u 2 ) , (61) 
FJl,u,3) = ^- - 4 + 8(10-n)w + O(w 2 ) . (62) 

Note that the first two terms of F±(l, u, 3) remain valid for general d < 4 due 
to the choice of the geometric factor A d in the definition of the renormalized 
coupling [14]. 
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6 Helicity modulus and superfluid density (stiffness constant) 

6.1 Bare theory 

The helicity modulus T measures the change of the free energy resulting from 
a small twist imposed on the order parameter [40,49]. This can be obtained 
from the change in the bare Gibbs free energy per unit volume F (r Q , uo, ho, k) 
[see Eq. (2)] caused by the k dependence of the field described by Eq. (15). 
For small k, this change is [40,49] 

AF = F (r , u , h , k) - F {r , u , h , 0) = \Tk 2 + 0{k 4 ) . (63) 

We are interested primarily in the physical case at h = for superfluid 4 He 
(n = 2). In the presence of a uniform superfluid velocity v s = h\t/m (relative 
to the motion of the normal fluid), phenomenological considerations within the 
two-fluid model identify the additional energy density as the kinetic energy 
density \p s v 2 = k B TAF . This implies that, for small k and h = 0, the 
superfluid density p s should be related to T by [40] 

Ps = k B T (m/hf T . (64) 

For 0(n) symmetric magnetic systems p s is the stiffness constant [38,60]. 

It turns out that at two-loop order our perturbative treatment is plagued by 
spurious Goldstone singularities for h = 0, T < T c . To circumvent this defect 
of perturbation theory we shall first work at h ^ and define the helicity 
modulus as 

(65) 

We shall see that the Goldstone singularities are cancelled in the limit h — > 0. 

In order to relate T to the Helmholtz free energy T we consider the Legendre 
transform r (ro, Uq, M , k) = Fo(r ,Uo,ho,k) + M ho- Differentiating F = 
Tq — Moho with respect to k 2 at fixed ho and making use of OTo/dMo = ho, 
we obtain (dFo/dk 2 )h = (dT /dk 2 )M and hence the bare helicity modulus 

(66) 

k=0 



d 

T = 2 ^ i7 o(n),wo, h , k) 
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where, in two-loop order, r is given by Eq. (26). Thus, from Eq. (66), we 
obtain 



T = Ml + + 2 Urn (^j + 0(4 ul In u„) , 

e=/^_ + /^_ "-I/ ^Bl 

J p r 0L +p 2 7 r 0T + p 2 k 2 7 (r 0L +P 2 )(^ t + P 2 ) 



(67) 
(68) 



where f 0L and f T are given by Eq. (24) and the two-loop diagrams repre- 
senting X are given in Fig. 3. The d = 3 poles enter only through the order 
parameter Mq since, being independent of k, they drop out of X upon dif- 
ferentiation. In d — 3, therefore, M 2 is given for ho — > by Eq. (41) while 
in O and X we may simply replace r by r , the difference affecting only 
terms of higher order. With this replacement, we write the one- loop contribu- 
tion as 0(ro, uo, M , k) = Q(r' ,uo, M , k) + 0(-Uq, Mq ln-u ); the corresponding 
two- loop contributions to the free energy are denoted by Xa, Xb, • • • , Xl (see 
Sec. 3 and App. C). 



X 



Both and 2dX /dk 2 \ k =o contain Goldstone singularities 
h — > below T c . Expanding up to O(u ) we obtain for r < 



h for 



0: 



.12tt 

@1 = - 77^7 M 



-2r{,) 1/2 + 0(a 1 / 2 ) 



+ ©i + O(m ,m In w ) , 



3a- 1/2 + n-12 + 0(a 1/2 ) 



= x: 1 (-2r + 2x: 1 )- 1 = ^(i-^)- 1 . 



(69) 

(70) 
(71) 



The other divergent terms come from diagrams C and I in Fig. 3 and are given 
(see App. C) by 



2d k 2X C 



U 



8tt 2 

2d k 2X! = -^u 



3 3 v 1 



(72) 
(73) 



These divergencies, however, are cancelled in the sum 

1 



lim 



0! + 2d k 2X c + 2d k2 X I 



16, n 

5 - n In 2 

3 



(74) 



The remaining terms are all finite for ho — > [see Appendix C] and their sum, 
along with Eq. (74), gives 
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T - i 2r '^ I —(-2r'Yl 2 - 1 (n I 2Wn L^IL^! 



1 



67T 2 



u [n — 15 — 3(n — 3) In 3] + 0(u Q , u lnu ) . 



(75) 



for ho — > and Tq < 0. Rewriting this in terms of the correlation length to 
absorb the logarithms of Uq we obtain finally 



8w £_ 24tt 
1 



2tt 2 



n — 3) In 3 



864tt 2 
urf- + 0(ug£) 



(1169n- 110) 



(76) 



for the helicity modulus in three dimensions. 



6.2 Amplitude function ofT in two-loop order 



The renormalization of the helicity modulus T(£_, uo, 3) or superfluid density 
Ps(£,-, u o,3) requires only the renormalization of the coupling constant Uq. In 
renormalized form, T R is obtained from T(£_,u ,d) according to 



u, /i, d) = T(£_, 2 A d 1 u, d) 

= £- d A d G(^_,u,d). 



(77) 
(78) 



The quantity of interest here is the dimensionless amplitude function G(l, u, d) 
which enters the final expression for T given below in Eq. (81). From Eqs. (76)- 
(78) with Z u (u, 1) and Z^(u, 1) in Eqs. (7) and (8) we obtain 



^-(2378 - 683n) + 8(n - 3) In 3 
54 



u + 0{u 2 ) (79) 



for the amplitude function in three dimensions. The first two terms agree with 
the previous [17,39] one-loop result. The solution of the RG equation 



[d^ + f3 u d u ] T R (£_, u,fi,d) = 0, 



(80) 



together with the choice = 1 for the flow parameter, yields the following 

(d = 3) representation for the temperature dependence of T near T c 



T = (4tt)- 1 C 1 G(1, «(/-), 3) 



(81) 



where the effective coupling u(l-) is determined by Eq. (A. 15). 
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6.3 Amplitude function of dx T (q) 1 /dq 2 in two-loop order 



For p s , there exists an alternative definition due to Josephson [41] in terms of 
the transverse susceptibility x T (q) at finite wavenumber q 



x t (q) 1 



g=0 



X T {q) = e^ x (0 o ,(x)0o,(O)) , 



(82) 
(83) 



where the system is considered to be in a homogeneous state with h ^ (that 
is, with k in Eqs. (15) and (18) set to zero). The inverse of x T (q) is equal to 
the transverse two-point vertex function 



r$(9) = r T + q 2 + 4m J G L (p) + 4(n + l)u J G T (p) 



- MulM* G L (\p + q\)G T (p) + ufo(q) + 0(4) 



(84) 



where u^(q) represents the negative sum of 22 two-loop digrams with two 
amputated transverse legs shown in Fig. 4. These diagrams are labelled by 
(1), (2), . . . , (22). The corresponding diagrams have also been given in Fig. 7 
of Bervillier [57] where, however, our diagram (8) is missing and where the 
diagram corresponding to our diagram (4) has an incorrect prefactor. [These 
errors do not affect dT t ^,/dq 2 \ q=0 , see Eq. (D.15).] Bervillier [57] calculated 
these diagrams within the e = 4 — d expansion. Here we shall evaluate them 
in three dimensions. 



We do not know of a rigorous proof of the equivalence between Josephson's 
definition of p S) Eq. (82), and the helicity modulus defined in Eqs. (63) and 
(64), but we shall verify this equivalence here in two-loop order by calculat- 
ing d~x T {q)~ l / dq 2 \ q= {) from Eq. (84). The result obtained for p s from Eq. (82) 
provides an important independent check of our calculation of the helicity 
modulus. Since the calculation is nontrivial we also provide some of the inter- 
mediate results in Appendix D. In particular, we present in Eqs. (D.6)-(D.ll) 
the Goldstone divergencies. They are more complicated than those of T but 
they are finally cancelled in the sum of the diagrammatic contributions. When 
expressed as a function of r' we denote T^^q) by T^\q). After a long calcula- 
tion, requiring considerably more computational effort than the calculation of 
T from the free energy, we obtain in three dimensions for h — > and r' < 
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+ (n - 1) + (n + 1) ^ + (n + 1) 2 - 

(1) (2) (3) (4) 



+ 256u 3 M$ 




+ 3 




+ (n + 1) 




+ (n-1) 




(5) 



(6) 



(7) 



(8) 



+ 512u 3 M% 



+ 3 




+ 2 




+ (n + 1) 




(9) 



(10) 



(11) 



(12) 



+ 128u 3 M* 






+ 2(n + l) + (n-1) 

(13) (14) (15) 



2048u 4 ,M o 4 




+ 2 




+ (n-1) 




(16) 



(17) 



(18) 



-4096m 4 M* 



- 32w 2 




(19) 



(20) 



(21) 



(22) 



Fig. 4. Two- loop diagrams for the two-point vertex function r o7 ;(g), Eq. (84), with 

longitudinal and transverse propagators Gl = and Gt = 1 — , Eq. (25). 

The diagrams (5), (12), (17), (18), (19) and (22) contain Goldstone singularities for 
h -»• [see Eqs. (D.6)-(D.ll)]. 



_d_ 

dq~ 2 



q=0 



1 + l- Uo (-2r' )- 1/2 + 18 - n - 24 In 3 )^(-2r )- 1 



(85) 



Eq. (85) contains no logarithms in uq since, at this order, there are no d = 3 
poles of dT^(q) / dq 2 \ q=0 . In terms of the correlation length £_ we obtain 
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d o 



r^ 2) (^,M ,g,3) 



g=0 



1 + ^-u £- + 8 - n - 12 In 3 



(86) 



Using Eqs. (82) and (42) this indeed reproduces T(£_, u , 3) as given in Eq. (76) 
and thereby proves Eq. (64) up to two-loop order. 

(2) ° (2) 

The renormalized vertex function 1^ is obtained from 1^ according to 

if } (£_, u, n, q, d) = Z+T%\i_, ^A^Z u Z7 2 u, q, d) . (87) 



On the basis of dimensional arguments we define the dimensionless amplitude 
function fx as [17] 



M^-,u,d) 



q=0 



Using Eqs. (87) and (88) and substituting Z u (u, 1) and Z^u, 1) we obtain in 
two-loop order in three dimensions 



f T (l,u, 3) = l + -« + 



48£ 



-An- 128 In 3 



u + 0(u 3 



(89) 



The first two terms agree with the earlier one-loop result [17,39]. Multiplying 
Eq. (89) by the amplitude function in Eq. (50) we indeed recover Eq. (79) 
from Eqs. (64), (81) and (82) according to 



G(l,u, 3) = 47r/ (l,«,3)/ T (l,«,3). 



(90) 



7 Results and discussion 



Within the minimally renormalized field theory in three dimensions we have 
derived the two-loop contributions to the amplitude functions of the follow- 
ing 0(n) symmetric quantities of the 0(n) symmetric 4 model at vanishing 
external field: 

(i) the square of the order parameter Mq = (0 O ) 2 , Eq. (50), 

(ii) the helicity modulus T, Eq. (79), 

(iii) the q 2 part of the inverse of the transverse susceptibility x T (q), Eq. (89), 
which yields the superfluid density p s (stiffness constant), Eq. (82), 

(iv) the specific heat above and below T c , Eqs. (61) and (62). 
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Goldstone singularities arising in an intermediate stage of the calculations have 
been shown to cancel among themselves. The resulting amplitude functions are 
applicable to the asymptotic critical region as well as to the non-asymptotic 
region well away from criticality (apart from corrections arising from finite- 
cutoff effects, from 6 terms and other higher-order couplings in H, Eq. (1), 
and from analytic terms). They provide the basis for 

(a) calculations of two- loop contributions to universal ratios Ai/A k and aija^ of 
leading and subleading amplitudes and a, which appear in the asymptotic 
representations [38] 

M = AJ^(l + a M |t| A + ...) , (91) 
Ps = A Ps \t\ (d - 2)u (l + a Ps \t\ A + ...) , (92) 

C ± = B + — \ t \- a (l + a ±\t\ A + . . .) , (93) 

(b) nonlinear RG analyses of non- asymptotic critical phenomena of 0(n) sym- 
metric systems above and below T c (for the general strategy of such analyses 
see Ref. [19]). 

These applications are of particular relevance to future experimental tests of 
the RG predictions of critical-point universality along the A line of 4 He (n = 2) 
[37]. 

The various amplitude functions are plotted in Figs. 5-9 vs the renormal- 
ized coupling u for the examples n = 1 and n — 2. In order to indicate the 
relative magnitude of the two-loop contributions we have also plotted the cor- 
responding zero- and one-loop approximations and, where available, the Borel 
resummed results (for n = 1 and n = 2 above T c [22] and for n = 1 below 
T c [23]). The curves terminate at the fixed points [16] u* = 0.0405 for n — 1 
and u* = 0.0362 for n = 2 (although extensions to u > u* may be needed in 
certain cases [21,44,45]). We comment on these curves as follows. 

For the superfluid density (n = 2), the one- and two-loop corrections for 
G(l,-u,3), see Fig. 5, each contribute about 10% of the zero-loop term l/8u* 
at the fixed point. (For n = 3 corresponding corrections are about 9% and 
5%, respectively, at u* = 0.0328.) The fact that the one- and two-loop contri- 
butions are of comparable magnitude suggests that higher-loop calculations 
including a Pade-type analysis or a Borel resummation are necessary before 
a reliable quantitative prediction can be made for amplitude ratios such as 
[17] Rj or a~/a Ps . A preliminary analysis [61] of experimental data [62,63] for 

p s and C ± , similar to the analyses in Refs. [19,39,64], indicate that the one- 
loop approximation for G(l,u*,3) is closer to the experimental result than 
the two-loop approximation. Since G(l,u, 3) can be considered as consisting 
of two factors Airf^ and f T according to Eq. (90), it is interesting to discuss 
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Fig. 5. Amplitude function G{u) = G(l,u, 3) for the superfiuid density (n = 2), 
Eq. (79), multiplied by u, as a function of the renormalized coupling u in zero-, one- 
and two-loop order. 

the latter amplitude functions separately. 



For the order parameter, we first consider for the case n — 1 for which the 
Borel summation result is known [23]. It is shown in Fig. 6 as the dot-dashed 
curve. As pointed out previously [23], the Borel result deviates only very little 
(by about 3% at the fixed point) from the zero- and one-loop result whereas 
the two-loop result is about 15% larger at u* '. Obviously, the leading order 
approximation is the better one in this case. We conjecture that this feature 
of the leading term will remain true also when n > 1 since the zero-loop term 
does not depend on n and since the one-loop term vanishes for general n [due 
to the choice of the geometric factor A d in Eq. (9)]. We note that for n = 2 the 
two-loop result for shown in Fig. 6 lies about 9% above the zero- and one- 
loop results, similar to the two-loop term for the case n — 1. Our experience 
with other amplitude functions above T c for which Borel summation results 
are available [22] is that the order in perturbation theory which is closest to 
the Borel results seems to be the same for n — 1, 2, 3. Clearly, in order to 
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n=1 




Borel 
(n=1) 



& 1 loop 
2 loop 







0.02 



0.04 



U 



Fig. 6. Amplitude function f^iu) = /</>(l, u, 3) for the square of the order parameter 
(n = 1, 2), Eq. (50), multiplied by u, as a function of the renormalized coupling 
u in zero-, one- and two-loop order. The dot-dashed curve is the Borel summation 
result given for n = 1 in Ref. [23]. 

make reliable quantitative predictions a higher-order calculation and Borel 
summation of the order parameter for n > 1 are necessary. 



Consider now the amplitude function fx of the transverse susceptibility shown 
in Fig. 7. For n — 2, the one-loop contribution at u* is about 10% whereas 
the two-loop correction is much smaller, being only about 1%. (For n = 3 
the corresponding corrections are about 9% and 1%, respectively.) We regard 
the smallness of the two-loop term as an indication of the quantitative reli- 
ability of the low-order formula for fx in Eq. (89). We infer from this and 
from our observations for that the considerable two-loop contribution to 
the amplitude function G is mostly due to the large two-loop term in and 
hence that the one-loop approximation for G is probably the most reliable at 
the present time. A substantiation of this conjecture by higher-order calcu- 
lations and Borel resummations would be highly desirable. From a practical 
point of view, the presumed reliability of the low-order result of fx is quite 
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Fig. 7. Amplitude function fr(u) = /t(1,«, 3) for dx T {q)~ l / dq 2 \ q= $ for n = 2, 
Eq. (89), as a function of the renormalized coupling u in zero-, one- and two-loop 
order. 

important since higher-order calculations of for n > 1 require considerably 
less computational effort than those of or of G. Thus, in a first step of 
future calculations of the amplitude function of the superffuid density, it will 
be sufficient to perform higher-order calculations only of f<f, for n > 1 before 
embarking on a long-term project of much more difficult higher-order calcu- 
lations of X T (q) or of T. We consider this important conclusion as a major 
result of our two-loop analysis. 

These considerations give some support to the good agreement obtained pre- 
viously [39] between the one-loop formula for the superffuid density and ex- 
perimental data [62,63] for 4 He in the (nonasymptotic) temperature range 
1CT 6 < t < 1CT 2 at several pressures near the A-line. A quantitative descrip- 
tion of the nonasymptotic region, however, depends not only on amplitude 
functions like G but also, crucially, on an accurate knowledge of the effective 
coupling u(l), which can be obtained from the experimentally determined spe- 
cific heat [14,64] . The formulas used to extract the nonuniversal initial value 
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Fig. 8. Amplitude function F_(u) = u, 3) for the specific heat below T c (n =1, 

2), Eq. (62), multiplied by u, as a function of the renormalized coupling u in zero-, 
one- and two-loop order. The dot-dashed curve is the Borel summation result given 
for n = 1 in Ref. [23]. 

of the effective coupling as a function of the pressure involve not only the RG 
exponent functions, which are known accurately from Borel summations [16], 
but also the amplitude functions of the specific heat F±(l,u, 3), which are 
given here to two-loop order in Eqs. (61) and (62). 

Figure 8 shows the amplitude function F_ of the specific heat below T c as a 
function of the renormalized coupling. As for the order parameter, we may use 
the known Borel summation results [23] for n = 1 (dot-dashed curve in Fig. 7) 
to infer the reliability of the low order approximations for n = 2. As noted 
previously [23], the one- loop approximation is the better one, differing from 
the Borel result at the fixed point by about 5% of the zero-loop term l/2u* 
compared with about 16% for the two-loop result. Evidently, the derivative of 
F_ at u — u*, which is needed for universal correction-to-scaling ratios [17], is 
also not well-approximated at two-loop order. Thus higher-order calculations 
of F_ for n > 1 are urgently needed. 
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Fig. 9. Amplitude function F + (u) = F + (l, u, 3) for the specific heat above T c (n =1, 
2), Eq. (61), as a function of the renormalized coupling u in one- and two-loop order. 
The dot-dashed curves are the Borel summation results given in Ref . [22] . 

In Fig. 9, we plot the amplitude function F + of the specific heat above T c , 
where Borel results are available for n > 1, as a function of the renormalized 
coupling [22]. Here, it is the two-loop rather than the one-loop approximation 
which is closer to the Borel results. This demonstrates that it is not clear a 
priori which (low) order of perturbation theory will provide the best approx- 
imation. 



Let us note finally that the Borel results [22,23] for F + (n > 1) and for F_ 
{n — 1) in three dimensions neglect the leading poles (in four dimensions) of 
the additive renormalization A(u, e) in Eq. (58) beyond two-loop order. Since 
a resummation of the (as yet unknown) higher-order terms in the RG function 
B(u) associated with A(u, e) (see Eqs. (2.32) and (2.33) in Ref. [17]) is expected 
to yield a small correction [65] of 0(77), we do not expect these Borel results 
for F± to be affected strongly by this approximation [22]. However, at the level 
of accuracy anticipated in future experiments [37], it is likely that the present 
uncertainties regarding B(u) will become significant, thus B{u) itself will also 



30 



be needed with the improved accuracy provided by a resummation of its high- 
order perturbation series. We recall that the function B(u) enters not only the 
formulas for the universal amplitude ratios [17,38] such as A + /A~ , Rj , a+/a~ 
and a~/a Ps but also the formulas needed to determine the effective coupling 
u(l) from the specific heat [14,19,64]. 

In conclusion, our new two-loop results for n > 1 below T c provide addi- 
tional motivation and specific information on the strategy and direction of 
higher-order calculations planned for future theoretical research parallel to 
the considerable effort on the experimental side to test the fundamental law 
of critical-point universality [37]. 



A Correlation lengths 



In this Appendix we derive Eqs. (13) and (14) and relate £± to t. Above T c , we 
obtain £ + from Eq. (12) and the two-point vertex function [16] T^\q, r , Uq) = 
X + (g) _1 . In two-loop order we have 



r^\q,r ,u ) =r + q 2 - -(n + 2)A d u Q r 1 e/z + ^(n + 2) 2 (1 - -)A 2 d u 2 r t 



l-e/2 



16 



|2„,2„l-e 



32(n + 2)-^% 
v ' (47r) d 



J (e) = J dx j dyjrj— 

o o 
i i 

J 2 (e) = dx dy 



Me) 
e-1 

ye/2-l 



Me)q 2 + 0(q 4 ) 



o o 
which leads to 



[1 - y + y(x - x 2 )] 2 -"/ 2 ' 
[1 - y + y(x - x 2 )}^' 2 



u y -* + o(4), (a.i) 

(A.2) 
(A.3) 



C 2 = r o {l - 4(n + 2)^jWo 6/2 + 16(n + 2) 
- a 2 {e)ulr Q e + 0{ul)} , 
a 2 (e) =32(n + 2)r(e)(47r)- d [j (e)(e - l)' 1 + J 2 (e) 



2e 



(A.4) 
(A.5) 



where J (l) = 2vr and J 2 (l) = 2tt/27. Inversion of Eq. (A.4) gives for £~ 2 ^ 



1 + 4(n + 2)-A d u C + + a 2 (e)(u C + ) 2 + O((u C + ) 3 )} , (A.6) 
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where 



77+2 

a 2 (6) = ^(e-l)- 1 [l + 0( e -l)] (A.7) 



determines the coefficient of the d = 3 pole of r 0c in Eq. (10). Subtracting 
^ r o — r o — r o [ see Eq. (11)] from Eq. (A. 6) and letting e — > 1 we obtain 



r = C 2 { 1 + ^(n + 2) Mo £+ + ^(n + 2)( Mo £+) 2 [1 + 541n(24 Mo £ + )] 



+ 



- -r(n + 2) ln24 -C(n) 



e=l 



7T 



K£+) S 



+ 0« In M ) 



(A.8) 



Here we have added and subtracted a logarithmic term 27r _2 (n+2)('Uo£+) 2 In 24 
in order to conform^ with the representation of r f in Refs. [11] and [23] for 
n = 1 in terms of the bare coupling g Q = 24w of Ref . [1 1] . Now a convenient 
choice for C(n) is such that the term in the square brackets of Eq. (A.8) 
vanishes, i. e. 



C(n) 



1 



(n + 2) 



7T 



Air 

l-C Euler + ln— -2 In 24 

y 



(A.9) 



This yields Eq. (13) and corresponds to the choice C\ = in Eq. (A2) of 
Ref. [11] forn = 1. 



Below T c , £_ is defined, according to Eqs. (3.1)-(3.6) of Ref. [17], as 
r'o = CX+KC, e) - ^£ 2 Z r (u{Q, e) - 5r , 



(A.10) 



(see also Eq. (A4) in Ref. [23]) where the function i2 +(wo£+, e) is defined above 
T c and represents the series of ro in integer powers of «o£+ whose first three 
terms are given in the square brackets of Eq. (A. 6). The effective coupling u(l-) 
as a function of the flow parameter Z_ is determined by Eq. (A. 15) below. For 
6 = 1, the function £l 2 i? +(wo£-, 1) ~~ ^ r o can be read off, up to two-loop order, 
from the right-hand side of Eq. (13) which represents £,+ 2 R 0+ (u t; + , 1) — 5r . 

1 See in particular Eq. (All) of Ref. [23]. Correspondingly Eq. (3.2) of Ref. [23] 
should read 5r = -3iio7r~ 2 (d - 3) _1 + C(l)u§ with C(l) given by our Eq. (A.9) 
for n = 1. In line 18 after Eq. (B5) of Ref. [23], " replacing -{A' d ) b " should read " 
replacing — (— A' d ) b ". 
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Thus it remains to rewrite Z r (u(l-), 1) as a function of wo£- (see Eq. (A7) of 
Ref. [23]). This leads to Eq. (14). 



The connection between £ ± and t is given by [16,17] 



at = Z- 2 Q+(l,u(l+),d) exp J Cr(«(0)y, * > 0> ( A - n ) 

f dZ 

-2a* = C 2 Q-(l,«(i-),d) ex P J Cr(«(0)j > t<0 > ( A - 12 ) 

i_ 

a = Z r (u, e) _1 a , (A. 13) 

Uu) = fid^nZ r (u,e)- 1 \ o , (A. 14) 

l±du(l ± )/dl ± =f3 u (u(l ± ),e), (A.15) 

where Q±(l, w, d) are the amplitude functions of the correlation lengths £±. 
They are related according to d) = 3 — 2Q + (l,u, d). A convenient 
representation for Q + (l,u, 3) is given in Eq. (3.5) of Ref. 



B Diagrams of Tq 1 ^ 



In this appendix, we show the cancellation of (spurious) Goldstone divergences 
when the order parameter is calculated from the one-point vertex function 
r^ro, Uq, M ) = ho as given by the sum of 1PI diagrams with one amputated 
external leg. These diagrams have been given to two-loop order by Bervillier 
in Fig. 5 of Ref. [57] and by Shpot in Eq. (17) of Ref. [51] where they were 
evaluated by use of an e = 4 — d expansion. Here we shall work at d = 3. The 
analytic expression, up to order 0(uq), is 



h 



j- = r + AupMl + 12m f G L {p) + 4(n - l)u J G T {p) + u 2 Y (B.l) 

p p 

where Fo = Yoa + Yob + ■ ■ ■ + Y 0i is the sum of the two-loop contributions 
Y 0a = -144 J G L ( Pl ) J G L (p 2 ) 2 , (B.2) 

Pl P2 

Y ob = -48(n -I) J G T (pi) J G L (p 2 ) 2 , (B.3) 

Pl P2 

y 0c = -16(n- 1) J G L ( Pl ) J G T (p 2 ) 2 , (B.4) 



Pl P2 
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Y od = -16(n 2 -I) J G T ( Pl ) J G T (p 2 ) 2 , (B.5) 

Pl P2 

Foe = -96 J J G L ( Pl )G L (p 2 )G L (\ Pl +p 2 \), (B.6) 

Pl P2 

Y of = -32(n - 1) J J G T ( Pl ) G T (p 2 ) G L (| Pl +p 2 |) , (B.7) 

Pl P2 

Y 0g = 3456u M* J jG L ( Pl ) 2 G L (p 2 )G L (\ Pl + p 2 \), (B.8) 

Pl P2 

Y 0h = 38A(n-l)u M 2 j J G T ( Pl ) G T (p 2 ) G L (| Pl + p 2 |) 2 , (B.9) 

Pl P2 

Y 0t = 25Q(n-l)u M 2 J J G T ( Pl ) 2 G T (p 2 ) G L (| Pl + p 2 |) (B.IO) 

Pl P2 

with G L (p) and G T ( P ) given by Eq. (25). Eqs. (B.1)-(B.9) agree with Refs. [51] 
and [57]. The prefactor 256(n — 1)u Mq in Eq. (B.IO) agrees with the prefactor 
of the corresponding diagram of Shpot [51] who corrected the corresponding 
prefactor of Bervillier [57]. 

The quantities Y 0e and F / contain d = 3 poles, but since Y 0e = 2(uIMq)~ 1 X 0H 
and Y of = 2(u 2 M 2 )- 1 (X 0I + X 0J ) for k = [see Eqs. (C.1)-(C3)], the sub- 
stitution r = r' + 5r in Eq. (B.l) leads to the cancellation of these poles in 
the same way as in Eq. (27). Thus, 

= 2u 2 f(r' ,u ,M ) (B.ll) 



lim 

r— >1 



5r + u 2 (Y e + Y f ) 



with f(r' ,UQ, Mq) given by Eq. (31). At this stage, one has the choice of 
performing the integrations before or after solving iteratively for Mq. Carrying 
out the integrations first, one is led to Eq. (34). Inverting first, one finds 



Ml = ^(-r + x T l ) + ^(-2r + ^f' 2 + ±(n - l)x T 1/2 + «o*i 

Uo[ - - ^ ^f(r' ,uo,x T ) + 0(u 2 ,ul]nu ) (B.12) 



Y~Y e -Y f 



where 



f(r' ,u ,x T ) = 



n 



+ 2, {-2r', + 3x T l Y 12 ,n-l, 1 + 2^ 



In 



71 A 



2AUr 



+ 



In 



7T Z 



, (B.13) 



with w given by Eq. (36), that is, / = / in Eq. (31) in lowest order in uq, and 
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n - l + 9w i/2 + 3w -i/2] | ( B14 ) 

is the 0(«o) contribution from the expansion of the integrals in Eq. (B.l) at 
one-loop order. In three dimensions we obtain for finite w 



Y a = 9(2n 2 )~\ 




(B.15) 


y 6 = 3(27r 2 )- 1 (n- l)w 1/2 , 




(B.16) 


Y c =(2n 2 )- 1 (n-l)w- 1 / 2 , 




(B.17) 


F d =(27r 2 )- 1 (n 2 -l) ! 




(B.18) 


Y 9 = 9(2n 2 )-\l-w), 




(B.19) 


Y h = 3(2n 2 r l (n-l)(l~2w l ^)(l- 


w)(l-4w) _1 , 


(B.20) 


% = n -*(n - 1)( W -V2 _ 2)(1 - «;)(! 


-41U) -1 . 


(B.21) 



Goldstone divergences ~ w^ 1 ^ 2 arise from Eqs. (B.14), (B.17), and (B.21) 
for ho — > 0. They cancel among themselves in Yi — |(Y C + Yj). Summing the 
remaining terms, one obtains the same result, Eq. (35), as derived via the free 
energy T. 



C Two-loop diagrams of r at k ^ 

The contributions X a, X b, ■ ■ ■ , X l in Eq. (26) at finite k are represented by 
the two-loop diagrams in Fig. 3. Diagrams A-G consist of products of exactly 
calculable one-loop integrals [24]. The integral expressions of the diagrams 
H L read 

X 0H = -48u 2 M 2 J J R L (p 1 )R L (p 2 )R L (\p 1 + p 2 \), (C.l) 

Pl P2 

X 0I = -lQu 2 M 2 J J R T ( Pl )R T (p 2 )R L (\ Pl + p 2 \), (C.2) 

Pl P2 

X 0J = -16(n-2)u 2 M 2 J jG T ( Pl )G T (p 2 )R L (\ Pl + p 2 \), (C.3) 

Pl P2 

Xok = -96u 2 M 2 J J R LT (pi) Rlt{P2) Rl(\Pi + p 2 |) , (C.4) 

Pl P2 

X 0L = -32u 2 M 2 J J R LT (pi) RtM Rt(\Pi + p 2 |) ■ (C.5) 

Pl P2 

Here we have used the elements of the matrix K^ 1 of k-dependent propagators 



Uo y 1 = J^{27 + (n-l) 
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Rl(p) = 
Rt(p) = 



1 R L (p) R lt (p)^ 

k Rtl(p) Rt(p) J 
ror +p 2 + k 2 

detAT(p,k) ' 
r 0L +p 2 + k 2 

detX(p,k) ' 



Rlt(p) 
Rtl{p) 



-2ikp 
det K(p, k) 

-Rlt(p) 



(C.6) 

(C.7) 
(C.8) 



where K(p,k), f L, f 0T , Gl[p) and G T (p) are given in Eqs. (23)-(25), respec- 
tively. The contributions dXj/dk 2 \ k=0 to the helicity modulus, Eq. (67), will 
be denoted by d k iXy The quantities d k iX c = d k 2[2u J Rl(Pi) f P2 Rt{P2)] 
and d k 2Xj are divergent for h — * 0. While d k 2X c is easily evaluated in terms 
of standard one-loop integrals, d k 2Xj is given, in three dimensions, by 



20*3 A". 



-4/3 



u 



(4tt) 3 
i l 

U = j dx j dy 



8(1 + a)h+ (l + 4a)/ 2 + 4a/ 3 - (6 + 8a)J 4 , (C.9) 



-1/2 



fi(x,y) 



[1 — y + y(# — x 2 )] 3 / 2 ' 



(CIO) 



/i(z,y) = m(y + a) - ln(l - y + a) , / 2 (x, y) = xy (xy + a) 1 , (C.ll) 
/ 3 (x,y) = lna-ln(l - y + a) , fa(x, y) = xy (1 - y + a)' 1 , (C.12) 

where a is given by Eq. (71). To illustrate the evaluation of these integrals, 
we consider J 4 . Thus, 



dx / dy 



, l l 2 x(l 



y + a) 



o o 

oo 1/2 

dz j du 
o 



[1 — y + y(x — x 2 )] 3 / 2 
(z - u 2 f/ 2 



1/4 

oo 



-It 



z- l/4 + a(* + 3/4)] 
(4z - l)- 1 / 2 



1/4 



^ [*-l/4 + a(z + 3/4)] 



2tt 



df 



(t 2 + l) {f + Aa{\ + a)- 1 ) 



3a 



-V2(l + a )l/2_ 2 



(C.13) 

(C.14) 

(C.15) 

(C.16) 
(C.17) 



In going from (C.13) to (C.17) we have used the substitutions z = y~ x — 3/4, 
u = x — 1/2, then -u = z 1 / 2 sin w and finally z — (t 2 + l)/4. While J4 exhibits a 
divergence for /i — * 0, the other contributions to d k 2Xj are finite in this limit 
and lead to Eq. (73). 



36 



The remaining contributions are finite for ho — > 0. For j = A, B and F they 
are readily evaluated in three dimensions as 



2d k 2X A = 2d k 2 3u (/ p R L (p)) = -(Stt 2 )- 1 ^ , 



3(87t 2 )- 1 m 



2<9 fc2 X B = 2«9 A;2 [3m (/ p i2 T (p) 
2,9 fc2 X F = 2,9,2 2(n - 2)« / Rt (Pi) J P2 Ct (P2 



(C.18) 
(C.19) 

7r 2 )- 1 (n-2)M . (C.20) 



9 fc2 X £ = <9 fc2 [2(n - 2)m /p iJ Rl(Pi)/ P2 G' t (p2)] vanishes for h -> 0, X G = 
n(n — 2)w [/ p Gr(p)] 2 is /c-independent, and X D = 4w [/ p -Rlt(p)] 2 = 0. For 
j = H, J, K and L we use splitting by partial fraction and obtain in three 
dimensions for ho — > 



29 fe2 X // = -(87r 2 )' 1 Mo , 
29 fc2 Xj = -(24vr 2 )- 1 (n-2) Mo , 

2 ( 9 fc2 X i , = 7r- 2 (l-ln3) Mo , 
2<9 fc2 X L = (37T 2 )^ 1 (21n2- 1)m - 



(C.21) 
(C.22) 
(C.23) 
(C.24) 



(2) /a„2 



D Contributions to dT Q rf/dq 



We denote the integral expressions of the two-loop contribution <&(q) in Eq. (84) 
by i — 1, 2, . . . , 22. The derivative of with respect to g 2 at q = yields 



d q ,F$(q) q=Q = l+ u V + u 2 d q ? <S>(q)\ q=0 + 0(4) , (D.l 

Grip) ~ G L (p) - (f 0L - r 0T )G L (p) 2 



8 _ er 0T + rfr 0L 
d (r L - for) 2 



(D.2) 



Expanding \I/ with respect to m an d replacing r by Tq we obtain in three 
dimensions 



u q = —u r \l + 0(a 1/2 )j + ufo 1 + O(u 3 ) , 

= -— 2 t 2 [Sa' 1 / 2 + n - 8 + O^ 1 / 2 ) 
^ = (-2r^ + 2x; 1 )- 1/2 . 



(D.3) 

(D.4) 
(D.5) 
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Eq. (D.4) exhibits a Goldstone singularity ~ a x l 2 with a given by Eq. (71). 
Further Goldstone singularities are contained in the two-loop diagrams (5), 
(12), (17), (18), (19) and (22) in Fig. 4. We find, in three dimensions, 



a,2$, = 7r V 



2_2 



a 



-1/2 _ 4 + 0(0,1/2) 



d q 2<& l2 = -2 (27vr 2 )- 1 (n + l)r 2 [cr 1 -91n(9a) - 12 + 0(a 1/2 ) 



^ 2 $ 17 = (6tt 



2\-l r 2 



12aT 1/2 - 7 - 80 In 2 + 0(a 1/2 ) 



18 



(27vr 2 )- 1 (n - l)r 2 la' 1 - 18 ln(9a) - 33 + 0(a 1/2 ) 



19 



(54tt 



2\-l 2 



4a- 1 - 72 ln(36a) - 87 + 0(« 1/2 ) 



( 9 <?2 $ 22 = (277r 2 )- 1 (n + l)r 2 a' 1 , 



(D.6) 
(D.7) 
(D.8) 
(D.9) 
(D.10) 
(D.11) 



where d q 2<&i = d$i/dq 2 \ q=0 . Summing up these contributions we get the finite 
result for h — > 



lim {^i + <9 ?2 [$ 5 + $i 2 + $i 7 + $ 18 + $i9 + $ 22 ]} 

ho— >0 

= (37r 2 )- 1 (-2r[ ) )- 1 [10 -4n- 48 In 2] . 



(D.12) 



The g-independent diagrams (l)-(4) and (13)— (15) do not contribute. The 
remaining diagrams give for ho — > in three dimensions 



a g2 $ 6 = 37r- 2 (-2r[ ) )- 1 , (D.13) 

^ 2 $ 7 = (n + l)7T- 2 (-2r;)- 1 , (D.14) 

<9 9 2$ 8 = 0, (D.15) 

^2$ 9 = -8 (3 7 r 2 )- 1 (-2r[ ) )- 1 , (D.16) 

^2$ 10 = -(27r 2 )- 1 (-2r[ ) )- 1 [l + 81n(3/2)] , (D.17) 

^2$ 11 = 7r- 2 (-2r;)- 1 [3-81n2] , (D.18) 

9 9 2$ 16 = (2 7 r)- 2 (-2r[ ) )- 1 [48 ln(3/2) - 1] , (D.19) 

<9,2$2 = -ii- 2 {-2r' Q )- 1 [1 + 16 ln(3/4)] , (D.20) 

( 9 9 2$ 21 = (127r 2 )- 1 (-2r;)- 1 . (D.21) 



This leads to Eq. (85). 
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